ChinaXiv 合 作 期 刊 


基于 概率 霍 夫 变 换 的 太阳 射 aili 1 型 暴 自动 识别 及 参数 提取 
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摘 要 : HEDERO E RAR TOA (Coronal Mass Ejection, CME) 期 间 产 生 的 非 热 高 能 粒子 

流 和 激 波 是 产生 灾害 性 空间 天 气 事 件 的 两 种 主要 因素 ， 它 们 主要 观测 特征 分 别 是 具有 快速 频率 
漂移 特征 的 射电 型 暴 和 具有 较为 缓慢 AIH YH DUR. AR EBD TARE A 
X # (Probabilistic Hough Transform 的 方法 ， 在 观测 数据 中 自动 识别 出 太阳 射电 了 H 型 暴 和 
了 型 射电 暴 ， 并 提取 射电 三 和 本 型 暴 的 特征 和 参数， 最 后 利用 识别 出 的 射电 本 型 暴 的 物理 参数 估 
计 了 CME 激 波 的 速度 。 
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作为 一 个 非常 活跃 的 天 体 ， 太 阳 上 经 常会 有 剧烈 活动 或 爆发 现象 〈 例 如 耀 斑 和 CME) o ME 
斑 和 CME 过 程 中 产生 的 高 能 粒子 和 激 波 等 在 日 地 空 间 中 传播 ， 有 可 能 会 到 达 地 球 ， 由 此 引发 一 
系列 的 日 地 物理 效应 (例如 太阳 高 能 粒子 事件 、 地 磁 暴 等 ) ， 对 人 类 的 航天 工程 、 通 讯 设施 、 
以 及 长 距离 输油管 道 等 都 会 有 非常 重要 的 影响 。 因 此 ， 为 了 对 空 x 间 天 气 进行 有 效 的 预报 和 预警 
减少 灾害 性 空间 天 气 对 人 类 的 生产 生活 的 破坏 ， 继 给 对 炮 兽 和 CWE 过 程 中 产生 的 高 能 粒子 和 激 
波 进行 深入 和 细致 的 研究 是 有 必要 的 。 

一 般 来 说 ， 天 体 射 电信 号 都 十 分 微弱 ， 但 太阳 射电 暴 的 信 号 略 强 。 太 阳 射 电 爆 发 指 的 是 
太阳 活动 期 间 在 射电 波段 出 现 的 剧烈 且 短 促 的 流量 增强 现象 ， 根据 其 在 射电 动态 频谱 上 展现 的 
不 同 的 形态 ， 可 以 分 为 太阳 射电 I-V 型 射电 暴 。 太 阳 射 电 I 型 和 开 型 暴 是 目前 研究 的 比较 多 的 
两 种 类 型 的 太阳 射电 暴 ， 分 别 是 由 泡 王 和 CME 期 间 产生 的 高 能 电子 束 流 和 激 波 所 导致 的 。 
具体 来 说 ， 泗 斑 和 CME 期 间 的 磁 重 联 加 速 电 子 形成 高 能 电子 ， 高 能 电子 束 流 以 亚 相 对 论 速 
FT 度 〈 约 为 0.2~0.6 c) vere Wk Be Ped E Jd FL Ss de aE ED, K 

. I 型 暴 在 射电 动态 频谱 上 展现 出 快速 的 频率 漂移 ， 射 电 I[ 型 暴 是 高 能 电子 束 流 的 最 佳 示 踩 器 。 
而 CME 在 日 晃 和 行星 际 空间 快速 运动 时 ， 当 CME 的 速度 超过 本 地 的 阿尔 芬 速 度 时 ， 会 产 
生日 晃 激 波 或 行星 际 激 波 。 目 前 通常 认为 太阳 射电 开 型 暴 是 由 于 激 波 在 日 冕 中 向 外 传播 引起 等 
离子 体 振荡 产生 朗 绪 波 ， 再 转换 为 电磁 波 辐射 并 以 本 征 等 离子 体 频率 和 二 倍 频 向 外 辐射 。 所 以 ， 
射电 开 型 暴 的 观测 特征 一 般 具 有 基 频 和 二 次 谐 频 结构 ， 由 于 激 波 的 运动 速度 远 低 于 高 能 电子 束 
流 ， 因 此 太阳 射电 开 型 暴 在 射电 动态 频谱 上 通常 表现 为 缓慢 的 频率 漂移 。 射 电 开 型 暴 可 以 被 视 
作 CME 激 波 的 最 佳 示 蹊 器 。 
目前 通常 认为 太阳 射电 英 型 和 开 型 暴 的 主要 辐射 机 制 为 等 离子 体 辐 射 ， 根 据 等 离子 体 辐射 
的 特点 可 以 通过 射电 暴 的 信号 频率 直接 得 到 射电 辐射 源 所 在 区 域 的 密度 。 再 根据 日 冕 密度 模型 
可 以 获得 辐射 源 所 处 的 日 冕 高 度 ， 进 而 获得 辐射 源 的 运动 速度 ， 例 如 CME 激 波 速度 ”。 
现 有 的 太阳 射电 观测 设备 普 忆 有 着 较 高 的 频率 分 辨 率 和 时 间 分 辨 率 ， 数 据 量 较 大 ， 如 果 仅 
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采用 人 工 识别 的 方法 的 话 ， 效 率 上 无 法 满足 实际 需求 。 到 目前 为 止 , 世界 上 已 经 发 展 了 多 种 太 
阳 射电 暴 识别 算法 。 这 些 算法 从 不 同 角 度 处 理 观测 数据 ， 如 对 流量 数据 使 用 阀 值 法 来 确认 爆发 
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开 型 暴 的 起 止 时 间 、 起 止 频 率 、 频 率 漂 移 率 等 重要 参数 ， 提 高 了 识别 的 准确 率 和 识别 速度 ， 并 
利用 识别 出 的 射电 暴 的 物理 参数 结合 相应 的 物理 模型 估计 了 CME 激 波 的 速度 。 这 不 但 对 空间 天 


气 的 预报 和 预警 有 重要 的 实用 价值 ， 并 对 研究 高 能 粒子 束 流 和 CME 在 日 揭 和 行星 际 空间 传播 物 


理 过 程 也 具有 重要 的 意义 。 


1 射电 I 型 暴 的 自动 识别 


程 ， 是 兆 席 /CME 等 太阳 爆发 事件 的 先兆 。 


AT Tere 
图 1 Phoenix 4 射电 频谱 仪 、Nancay 10 米 波 射电 频谱 仪 和 
WIND/WAVES 射电 探测 器 在 2014 年 2 月 16 日 获得 的 一 例 
超 宽 带 射电 型 暴 的 综合 频谱 图 


射电 型 暴 是 一 种 最 常见 的 太阳 射电 欢 ， 是 太阳 剧烈 活动 后 产生 的 非 热 融 能 电子 束 在 日 园 
和 行星 际 空间 中 运动 的 最 佳 示 踪 器 ， 同 时 也 能 最 先 反 映 太 阳 暴 能 量 释 放 和 粒子 加 速 等 物理 过 


Figure. 1 A synoptic spectrogram of a type II] radio burst with super-wide 
frequency band observed by Phoenix 4 decimeter spectrometer, Nancay 


decameter array and WIND/WAVES on 2014 February 16°” 
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此 ，10MHz 以 下 的 射电 II 


型 暴 的 频率 漂移 速度 要 远 低 于 10MHz 以 上 的 。 而 高 于 10MHz 的 射电 I 贡 型 暴 具 有 类 似 直 线 的 结构 ， 
这 也 使 得 本 文 使 用 概率 霍 夫 变换 算法 来 进行 射电 I 贡 型 暴 识别 成 为 可 能 。 


本 文 使 用 法 国 南 希 的 十 米 波 天 线 阵 (NDA, Nancay Decameter Array) 的 太阳 射电 观测 数 
据 进 行 II 型 暴 识别 ， 其 时 间 分 辨 率 为 1s， 其 观测 频率 范围 为 10-80MHz， 其 中 均匀 分 布 着 400 


个 通道 。 
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为 了 识别 出 观测 数据 中 所 包含 的 太阳 射电 暴 的 数据 ， 需 要 对 观测 数据 进行 处 理 。 处 理 过程 
分 为 三 个 部 分 : 

(1) 原始 动态 频谱 图 的 分 段 及 背景 去 除 ，(2) 动态 频谱 图 的 前 景 提取 及 噪声 处 理 ，(3) 识别 
器 型 暴 事件 并 提取 爆发 相关 信息 。 
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图 2 MERRIE 
Figure.2 Identification process of type II radio burst 


1.1 太阳 射电 频谱 数据 预 处 理 

在 太阳 射电 动态 频谱 图 上 ， 了 UI 型 暴 表 现 为 具有 快速 频率 漂移 的 信号 ， 且 不 同型 暴 事 件 的 
相对 强度 有 很 大 差异 ， 而 背景 信号 强度 则 在 几 小 时 的 时 间 内 逐渐 变化 。 因 此 我 们 把 原始 数据 切 
分 成 若干 小 片段 ， 这 样 有 助 于 去 除 背 景 效 应 ， 同 时 也 能 尽量 避免 信号 强 的 射电 暴 和 弱 的 射电 暴 
在 同一 个 数据 片段 中 而 导致 弱 的 射电 暴 不 被 识别 ， 经 过 多 次 实验 测试 ， 我 们 得 出 每 个 片段 为 
10 分 钟 的 观测 数据 可 以 有 助 于 射电 暴 的 识别 与 参数 提取 。 

在 NDA 可 观测 的 频率 范围 内 ， 我 们 将 每 个 数据 段 的 时 间 长 度 设 置 为 600s， 采 样 间 隔 500s, 
即 两 个 相 邻 数据 段 的 重 钱 时间 为 100s， 这 种 切 分 方式 可 以 保证 在 某 个 片段 中 恰好 被 切断 的 爆 
发 事件 在 下 一 个 片段 中 可 以 被 识别 。 因此 ， 切 分 之 后 每 一 个 数据 段 都 是 一 个 [600X250] 的 二 


维 数组 1(ti,f)。 
由 于 观测 设备 的 不 同 通道 之 间 存在 明显 的 增益 差异 ， 由 此 产生 的 通道 效应 在 频谱 图 上 表现 
为 一 条 条 高 亮 的 横 条 纹 ， 这 些 横 条 纹 会 造成 爆发 结构 间断 ， 有 时 甚至 会 这 盖 较 弱 的 爆发 事件 ， 


影响 对 太阳 射电 骏 的 识别 ， 所 以 必须 先 消除 观测 数据 中 的 通道 效应 。 
有 基体 处 理 过 程 为 ， 在 各 通道 的 数据 上 减 去 该 通道 所 有 数据 的 均值 ， 以 去 除 背 景 ， 只 保留 来 
自 太阳 的 射电 信号 ， 实 现 各 通道 数据 的 归 一 化 ， 有 具体 公式 如下: 
nal tw fj 
ref= hem -ZEA 


其 中 ，1(ti,) 与 1(ti, 友 ) 分 别 表 示 经 过 处 理 之 前 和 之 后 的 数据 ，a 与 5 分别 表示 该 数据 段 的 
开始 时 刻 与 结束 时 刻 的 索引 。 需 要 对 数据 段 I 中 的 每 个 值 都 进行 公式 1 所 示 的 计算 ， 以 生成 新 
的 消除 通道 效应 后 的 7'。 
1.2 Yen 算 法 二 值 化 及 图 像 腐蚀 
在 对 太阳 射电 动态 频谱 图 进行 预 处 理 后 ， 还 需要 对 图 像 进行 二 值 化 ， 将 图 像 灰 度 值 置 为 0 
或 255， 经 过 归 一 化 后 即将 灰 度 值 置 为 0 或 1。 使 其 目标 前 景 与 背景 被 区 分 开 来 ， 以 便 避 开 背 景 
的 影响 ， 利 于 寻找 目标 前 景 中 的 太阳 射电 三 型 暴 事 件 。 
可 以 对 太阳 射电 动态 频谱 图 进行 二 值 化 的 算法 有 很 多 种 ， 但 效果 各 有 优 劣 。 本 文 使 用 多 种 


自动 确定 阔 值 的 算法 对 观测 数据 进行 测试 ， 生 成 二 值 化 图 像 ， 原 始 图 像 与 使 用 各 种 效 值 法 进行 
二 值 化 的 图 像 如 图 3 所 示 。 
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图 3 Al A sy ce BS AY FEE eR 
Figure. 3 Binarization effect of various automatically determining threshold algorithm 


在 图 3 中 ，0riginal 代 表 原 始 图 像 ， 其 余 图 像 代表 各 种 算法 (如 : Isodata, Li, Mean, 
Minimum, Otsu, Triangle, Yen) 对 原始 图 像 进行 二 值 化 的 输出 结果 。 其 中 最 为 典型 的 大 津 法 
( 即 0tsu 法 ) ” ， 也 被 称 为 最 大 类 间 方 差 法 ， 就 是 一 种 使 用 较为 广泛 的 自 适 应 阔 值 化 方法 。 
算法 通过 不 断 尝 试 将 各 种 值 作为 闷 值 ， 来 计算 图 像 的 类 间 方 差 ， 选 取 类 间 方 差 最 大 时 的 冰 值 件 
为 对 该 图 像 进行 二 值 化 的 标准 。 此 时 ， 原 图 像 以 该 闷 值 作为 界限 ， 将 图 像 中 的 所 有 像素 分 为 背 
景 像素 和 目标 像素 两 大 类 
经 测试 ， 可 以 从 图 3 的 测试 结果 中 看 出 ， 只 有 Yen 算 法 能 较 好 地 确定 阐 值 ， 在 二 值 化 后 明确 
区 分 图 像 前 景 与 背景 ， 并 较 好 地 保留 了 前 景 的 信息 ， 即 爆发 的 信息 ， 未 与 背景 信息 混淆 在 一 起 。 
而 其 他 算法 的 输出 结果 未 能 准确 区 分 图 象 前 景 与 背景 ， 导致 爆发 信息 被 干扰 ， 难以 识别 。 

Yen 算 法 ” 是 一 种 自动 多 门限 闵 值 法 ， 能 够 根据 太阳 射电 频谱 图 像 本 身 的 灰 度 特性 ， 将 图 
像 分 割 为 目标 前 景 与 背景 两 个 部 分 。 与 0tsu 法 类 似 ，Yen 算 法 也 能 上 自动 选取 最 佳 闵 值 7， 实 现 图 
BREDA, ae Be EB ELT BY CH a Be AC EU, ES EE AE PEEL, OR EY 
景 和 目标 所 共同 提供 的 信息 量 最 大 。 
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图 4 经 过 二 值 化 与 iL 


Figure.4 Image after binarization and erosion 


WEAR, ERAR SY Hoh AS iS RET EO» ibe BET I BE TT SA RE, UE 
除 图 像 中 约 每 小 时 一 次 的 定 标 信号 ， 以 及 某 些 在 频谱 图 中 所 占 面 积 较 小 的 瞬时 干扰 信和 号。 


腐蚀 (erosion) 是 一 种 基本 的 形态 学 运算 , 可 以 寻找 图 像 中 的 极 小 区 域 ， 将 图 像 中 的 高 亮 区 
域 进 行 缩减 ， 在 其 输出 的 图 中 ， 高 亮 区 域 比 万 图 的 高 亮 区 域 更 小 。 
图 4 中 可 以 看 出 ， 利 用 腐蚀 算法 ， 频 谱 图 中 的 定 标 信号 与 瞬时 干扰 都 被 较 好 地 去 除 掉 了 。 
同时 ， 受 到 腐蚀 算法 的 影响 ， 爆 发 事件 的 结构 的 边缘 也 遭 到 腐蚀 ， 但 整体 结构 与 频率 漂移 率 变 
化 不 变 ， 不 影响 后 续 的 识别 。 
[3 BERLE ea ST ESE JAR BES LAR SE 
霍 夫 变换 是 一 种 特征 检测 方法 ， 被 广泛 地 应 用 在 图 像 分 析 、 计 算 机 视觉 等 领域 。 使 用 霍 夫 
变换 的 方法 ， 可 以 在 图 像 中 寻找 茶 些 图 形 ， 如 :; Hk Dl. ALB 
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i r 
transform principle in polar ordinate system 


图 5 ERANU” 


Figure.5 Hough transform principle™ 
常用 的 霍 夫 变 换 检测 直线 的 方法 , 是 运用 下 式 在 图 像 空 间 和 参数 空间 之 间 , 将 图 像 中 的 曲线 
转换 为 参数 空间 的 一 点 , 建立 对 偶 变 换 。 
p=xcosa+ysina #(2) 
AF, p 为 极 径 ，a 为 极 角 , aHLO~180° ; x 为 像素 点 相对 图 像 原点 的 行 坐标 ，y 为 像素 点 相 
对 图 像 原 点 的 列 坐 标 。 


图 6 霍 夫 变换 的 效果 ” 


Figure.6 The effect of Hough transform 

如 图 6 所 示 ， 输 入 的 图 像 被 从 图 像 空间 转换 到 参数 空间 后 ， 在 参数 空间 中 存在 明显 的 三 个 

亮点 ， 而 参数 空间 的 亮点 则 代表 着 输入 的 图 像 中 的 直线 结构 。 利 用 霍 夫 变换 的 直线 检测 功能 
我 们 可 以 在 太阳 射电 频谱 图 像 中 寻找 拥有 类 似 直 线 的 结构 的 太阳 射电 暴 。 

本 文 使 用 的 概率 霍 夫 变 换算 法 利用 解析 曲线 的 几何 特性 ， 尺 量 避 免 对 图 像 中 的 每 个 非 零 像 


[9] 


素 进行 计算 量 大 的 投票 过 程 ， 而 是 从 点 集中 随机 选取 一 个 像素 点 进行 分 析 。 

概率 霍 夫 变换 算法 具体 步骤 如 下 : 

(1) 每 个 区 间 对 应 一 个 累加 器 4( 9)， 从 每 个 区 间 的 点 集中 随机 地 选取 一 个 像素 点 ， 在 参 
数 空间 中 计算 9 对 应 的 r 值 ， 并 在 对 应 的 累计 器 A(r,0)_ANH1; 

(2) 从 点 集中 删除 该 点 并 更 新 累加 器 ; 

(3) 若 累 加 器 的 值 大 于 提前 设 定好 的 阐 值 ， 则 表示 已 检测 到 直线 ， 然 后 删除 该 直线 上 所 有 
的 点 ; 

(4) 重复 以 上 3 个 步骤 ， 直 至 点 集 为 空 

根据 太阳 射电 三 型 暴 的 频谱 信号 特征 ， 可 以 知道 爆发 事件 在 太阳 射电 动态 频谱 图 中 近似 一 
条 倾斜 的 直线 。 直 线 的 倾斜 程度 与 太阳 射电 I 型 暴 的 频率 漂移 率 相关 。 

对 太阳 射电 动态 频谱 图 使 用 霍 夫 变换 来 寻找 [型 暴 。 如 果 在 动态 频谱 数据 片段 中 存在 符合 
约束 条 件 的 直线 ， 则 可 以 认为 该 动态 频谱 数据 片段 中 可 能 存在 II 型 暴 事件 。 但 是 ， 在 数据 段 
可 能 还 存在 其 他 类 型 的 射电 暴 ， 例 如 太阳 射电 开 型 暴 ， 其 在 频谱 上 也 可 能 具有 类 似 直 线 的 图 像 
特征 。 所 以 需要 根据 识别 出 的 频率 漂移 率 来 限制 及 过 滤 其 他 类 型 爆发 事件 。 

动态 频谱 数据 片段 经 过 前 文 所 述 的 二 值 化 、 腐 蚀 处 理 后 ， 对 于 每 一 帧 二 值 化 后 的 图 像 进 行 
霍 夫 变换 ， 检 测 图 像 中 的 直线 结构 。 同 时 ， 已 知 在 NDA 观 测 数据 的 频率 范围 内 ， 太 阳 射 电 II 型 
暴 的 频率 漂移 率 约 为 1-10MHz/s， 所 以 需要 对 检测 出 来 的 直线 进行 筛选 ， 只 留 下 频率 漂移 率 符 


合 I 型 暴 的 直线 。 
若 经 过 筛选 后 ， 数 据 片段 
件 。 


同时 还 可 以 获取 线段 两 端的 多 


中 仍 存在 一 


图 7 霍 夫 变换 算法 识别 出 的 爆发 事 


Figure. 8 Burst event identified 
如 图 7 所 示 ， 图 中 有 直线 特征 的 结构 被 识别 
“ 标 值 ， 即 获取 爆发 
频率 。 图 8 中 所 示 事 件 的 起 始 频 率 为 50MHz， 终 1 


条 或 多 条 直线 ， 则 和 有 


件 


by Hough transform algorit 


和 定 该 数据 片段 中 包含 至 少 一 个 事 


nm 


出 来 并 标 红 ， 此 处 很 可 
事件 的 起 始 时 间 、 


率 为 -1. 5MHz/s， 负 值 表 示 II[ 型 暴 从 高 频 向 低频 漂移 。 
1.4 NDA 数 据 中 II 型 暴 的 识别 结果 


起 始 
上 频率 为 27MHz， 持 续 时 间 为 15 秒 ， 


能 是 一 个 型 暴 事 件 。 
频率 、 终 止 时 间 、 终 目 
频率 漂移 速 


使 用 以 上 识别 算法 可 以 高 效 地 在 10MHz 以 上 的 频谱 数据 中 识别 形态 类 似 直 线 的 本 型 暴 事件 ， 
并 获得 关于 其 起 止 频率 、 起 止 时 间 ， 频 率 漂 移 率 等 信息 。 为 之 后 对 爆发 事件 的 研究 提供 支持 。 


=> 


图 8 识别 出 的 不 同 参数 的 爆发 事件 
Figure.8 Burst events with different parameters identified 
如 图 8 所 示 ， 图 8 (a) 中 爆发 事件 的 起 始 频率 为 67MHz， 终 止 频率 为 23MHz， 持 续 时 间 为 17s， 
频率 漂移 率 为 -2. 6MHz/s。 图 8(b) 中 爆发 事件 的 起 始 频率 为 37. 3MHz， 停 止 频率 为 28MHz， 持 续 
时 间 为 6s， 频 率 漂移 率 为 1. 55MHz/s。 以 上 两 例 爆发 事件 的 识别 结果 证 明 识 别 算法 可 以 识别 从 


各 频率 开始 的 爆发 事件 。 

图 8(c) 中 爆发 事件 的 起 始 频率 为 58MHz， 终 止 频 率 为 46MHz， 持 续 时 间 为 6s， 频 率 漂 移 率 为 
-2. 1MHz/s。 图 8(d) 中 爆发 事件 的 起 始 频 率 为 48MHz， 终 止 频率 为 26MHz， 持 续 时 间 为 18s， 频 率 
漂移 率 为 -1. 22MHz/s。 以 上 两 例 爆发 事件 的 识别 结果 证 明 识 别 算法 可 以 识别 各 种 长 度 的 持续 时 
间 的 爆发 事件 。 


Freq/MHz 


09:45:39 09:46:29 09:47:19 09:48:09 09:48:59 09:49:49 


Time/UT 


图 9 识别 出 的 低频 的 爆发 事件 
Figure.9 Burst event with low frequency identified 
但 由 于 NDA 观 测 的 频率 范围 有 限 ， 且 在 观测 范围 内 的 低频 频段 中 存在 较 强 的 持续 不 断 的 干 
扰 ， 如 果 爆 发 事件 的 终止 频率 较 低 ， 这 些 干扰 信号 可 能 会 影响 算法 对 爆发 事件 终止 频率 的 识别 。 
如 图 9 所 示 ， 算 法 识别 出 的 爆发 事件 的 终止 频率 与 干扰 信号 所 在 频率 紧密 结合 ， 这 样 在 识别 爆 
发 事件 终止 频率 时 ， 会 有 一 定 的 误差 。 


2 射电 开 型 暴 的 自动 识别 


从 太阳 射电 观测 数据 中 识别 太阳 射电 开 型 暴 的 方法 与 识别 太阳 射电 匡 型 暴 的 方法 类 似 ， 需 
要 对 数据 进行 图 2 所 示 的 处 理 过 程 。 但 在 使 用 霍 夫 变 换算 法 识别 出 频谱 图 中 的 直线 结构 及 其 起 
止 点 坐标 后 ， 需 要 根据 开 型 暴 的 物理 特征 进行 筛选 。 例 如 :， 开 型 暴 的 基 波 所 在 的 起 始 频率 通常 
不 会 很 高 〈 通 常 <300MHz) ， 其 频率 漂移 速率 较为 缓慢 〈 通 常 <1 MHz/s) o 

太阳 射电 开 型 暴 的 主要 辐射 机 制 为 等 离子 体 辐射 ， 所 以 可 以 通过 射电 暴 的 信和 号 频率 直接 得 
到 射电 辐射 源 所 在 区 域 的 密度 。 根 据 等 离子 体 辐射 的 特点 ， 观 测 频率 fobps、 等 离子 体 辐射 频率 
fpe 和 电子 密度 ne 的 关系 可 表示 为 fobs = Sfpe> fpe(kHz) = 8.98Vne(cm-3)。 这 里 s 代 表 谐 波 数 ， 
s = 1 时 代表 基 频 结构 (Fundamental) ，s = 2 代表 是 二 次 谐 频 结构 (Harmonic) 。 在 此 基础 之 
上 ， 还 可 以 根据 日 冕 密度 模型 计算 射电 辐射 源 所 在 的 日 冕 高 度 或 位 置 ， 并 进一步 估算 出 CME 


nee ` Al AW He » $ sH [10] 
HEE. AEA Se ay BE TE EP BAI H KHAIRA — iE Newkirk (1961) 模型 
Ne = No X 10432RsAR #(3) 
其 中 , No = 4.2 x 10 4cm 性; Rs 是 太阳 半径 。 


基于 以 上 自动 提取 算法 和 日 冕 密度 模型 (公式 3)， 我 们 对 于 云南 天 文 台 米 波 太阳 射电 频谱 
仪 2012 年 11 月 13 日 的 观测 数据 进行 处 理 ， 识 别 其 中 的 开 型 暴 ， 并 提取 起 止 频率 、 起 止 时 间 、 


频率 漂移 率 等 参数 ， 并 计算 所 对 应 的 CME 的 特征 参数 (如 图 10)。 


Freq/MHz 


CME 特征 参数 提取 
识别 目标 ， 开 型 射电 暴 CH) * 
F 始 时 间 : 02:02:20 UT 

结束 时 间 : 02:05:44 UT 

开始 频率 :388 MHz 
结束 频率 : 158 MHz 

频率 漂移 率 : -0.56 MHz/s 
CME 激 波 速度 : 870km/s 

注 : * 中 的 H (Harmonic) 代表 开 型 射 
电 暴 的 二 次 谐 波 结构 

02:03 ooe * * 中 的 负 号 代表 开 型 射电 暴 的 频率 漂 
iii 移 是 从 高 频 到 低频 

图 10 利用 观测 数据 提取 CME 特 征 参数 


M4 


= 


Figure. 10 Extracting CME feature parameters from observation data 


TELE ERA RATE a 
02:05:44 UT 时 刻 的 速度 为 870k 


该 CME 的 速度 为 851km/s 数据 来 源 : 


https://cdaw. gsfc. nasa. gov/ 


该 型 射电 暴 的 物理 参数 ， 并 得 出 对 应 的 CME 在 02:02:20 UT- 
m/s, Tifa, LASCO C2A #1 7£02:24:06 UT 发 现 该 CME， 并 获得 


CME list/UNIVERSAL/2012 11/univ2012 11.html ) ， 通 过 对 比 


发 现 ， 这 两 种 方法 所 给 出 的 CME 
识别 出 的 参数 是 合理 可 信 的 。 


3 总 结 与 展望 


速度 非常 接近 ， 因 此 ， 我 们 认为 利用 概率 霍 夫 变换 算法 所 自动 


本 文 提出 了 一 种 基于 概率 霍 夫 变换 算法 的 太阳 射电 II 型 和 开 型 暴 识 别 算法 。 首 先 通过 一 系 


图 像 中 的 直线 结构 ， 获 取 IIL 和 I 


这 些 信息 进一步 得 选 频 率 漂 移 率 符合 IJ 型 和 开 型 暴 特 征 


列 的 算法 对 太阳 频谱 图 像 进行 预 处 理 : 多 次 实验 测试 选择 合适 的 时 间 长 度 对 原始 观测 数据 进行 
合理 切 分 ， 对 各 通道 进行 均值 消去 ， 消 除 通道 不 均匀 性 的 干扰 ， 对 图 像 进行 二 值 化 ， 提 取 前 景 
信息 ;对 图 像 进行 腐蚀 ， 消 除 一 些 窗 带 的 射频 干扰 信号 。 然 后 使 用 概率 霍 夫 变换 算法 寻找 频谱 


I 型 暴 事件 的 起 止 时 间 、 起 止 频率 和 频率 漂移 率 等 信息 ， 并 根据 
直线 结构 。 


本 文 还 利用 该 方法 从 云南 天 文 台 2012 年 11 月 13 日 工 型 暴 的 观测 数据 中 自动 提取 了 [型 暴 的 


物理 参数 ， 并 结合 相应 的 日 蝎 密 度 模型 获得 CME 的 特征 


T 


参数 。 这 些 参数 对 空间 天 气 的 预报 预 


警 具有 实用 价值 。 该 方法 虽然 可 以 识别 典型 的 开 型 又 ， 但 由 于 开 型 暴 容易 与 其 他 类 型 射电 暴 同 


好 的 、 对 多 种 数据 来 源 普遍 适 月 


时 发 生 ， 信 号 有 混 刁 ， 未 来 我 们 将 在 更 多 的 工 型 暴 数据 上 对 该 算法 进行 检验 和 改进 ， 以 达到 更 


的 识别 效果 。 


实验 结果 表明 ， 该 方法 能 有 效 地 识别 太阳 射电 匡 型 暴 与 工 型 暴 ， 并 获取 射电 姑 的 起 止 频率 、 


起 止 时间 和 频率 漂移 率 等 信息 ， 
同时 ， 概 率 霍 夫 变 换算 法 上 


并 能 快速 地 对 长 期 、 大 量 的 观测 数据 进行 识别 和 分 析 。 
标准 霍 夫 变换 算法 的 计算 量 更 小 ， 计 算 过 程 更 短 ， 更 适合 在 观 


测 过 程 中 实时 识别 太阳 射电 I 型 和 开 型 暴 ， 提 取 观 测 参数 得 出 CME 激 波 速度 等 关键 物理 信息 ， 


实现 空间 天 气 预报 预警 ， 维 护 本 


日 关 人 员 、 设 备 的 安全 。 
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Abstract:Nonthermal energetic particle flow and shock wave generated during flares and coronal 
mass ejection are the two main factors that produce disastrous space weather events. Their main 
observational characteristics are type III radio burst with fast frequency drift feature and type II radio 
burst with relatively slow frequency drift feature. 

This paper mainly introduces the method of using the probabilistic Hough transform to 
automatically identify the type III and type II solar radio burst in the observation data, and to extract 
the characteristic parameters of the type III and type II radio burst, finally the velocity of the CME 
shock wave is estimated by using the physical parameters of the identified radio type II bursts. 
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